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Abstract 

We consider the asymptotic behavior of the (one dimensional) two-species annihilation reaction 
A + B — !■ 0, where both species have a uniform drift in the same direction and like species have a 
^—^ I hard core exclusion. Extensive numerical simulations show that starting with an initially random 

^— ^ ■ distribution of A's and B's at equal concentration the density decays like t^^/'^ for long times. 
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This process is thus in a different universality class from the cases without drift or with drift in 
different directions for the different species. 
PACS numbers: 5.40.+J, 82.20.Mj, 2.50.-r 
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Ph The irreversible two species annihilation reaction A + B ^ has been studied for quite some 

o' 

time as an example of a reaction diffusion process where fluctuations are important, so that the 
^ , density decays more slowly than one would predict from the mean field rate equations (at least in 

low dimension) ||l], |2|, |3|, ^, |5|, |6| . The standard picture |j|] is that in a region of size L it takes a time 
on the order of L^ for all the particles to react, since they must diffuse around the region in order to 
annihilate. The remaining density will be proportional to the intial excess of either type A or type 
B particles in this region, which is proportional to the square root of the volume. Thus at time t one 
expects that the concentration c{t) will behave as 
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it) ~ [c(0)L'^] ' /l'^ ~ [c{0)t''/^] ' /t'^/2 = c{Qflh-'"\ (1) 
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for dimension d < 4, i.e. c(0)^'^t ^'^ in dimension one, which has been verified rigorously |^. In 
comparison the mean field result is c(t) = [t + l/c{0)]~^. 

In order for this picture to be valid, it is necessary that the distribution of particles be actually 
determined by diffusion. The addition of a drift field can invalidate this assumption. 

That this happens when the two species drift in opposite directions is not unexpected. If the 
two species have a relative velocity v, then using the same reasoning as above, one determines the 
concentration to be 0, ^ 

c(t) ~ [c(0)z;t(^+i)/2] '^' / [^t('^+i)/2] = [c(0)^]i/2r ('^+1)/^ (2) 

for dimension d < 3, or [c(0)u]^'^t^^'^ in one dimension. 

What is perhaps unexpected is that if the two species have a drift in the same direction, the 
result may still be different from the expected t~ ''^. This is because even after subtracting off the 
average motion, one does not necessarily recover diffusive behavior. 

In order to proceed, we must specify the model further. To keep matters simple we consider 
a strong drift — particles move either to the right or not at all. We consider two types of particles 
moving on the one dimensional lattice (with periodic boundary conditions). At each (micro-)time 
step, we randomly pick a site. If that site is occupied, we attempt to move the particle to the right. 
If the site on the right is unoccupied, the jump succeeds. If the site is occupied by a particle of 
the same species, the jump fails. If the site is occupied by a particle of the opposite species, both 
particles annihilate. 

We now see why the diffusive picture might be invalid. If we just considered one species of 
particles, the model described above is the Asymmetric Simple Exclusion Process (ASEP), and the 
long time behavior is not governed by a linear diffusion equation. Instead one must consider a 
nonlinear stochastic equation such as the noisy Burgers equation p, 

^ + n^ = u^ + — (3) 

dt dx dx"^ dx^ 
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Figure 1: Log-log plot of density vs. time for system of size 4 x 10 . Initial densities are 0.9, 0.5 and 
0.2. 



where p is a rescaled density, v is the viscosity (representing the lattice spacing) and ^ is a random 
noise term, e.g. white noise where the covariance is 



{i{x, t)i{x', t')) = 2u5{x - x')5{t - t'). 



(4) 



Thus the exclusion rule, even for very low densities, can play a crucial role in the dynamics of the 
annihilation process, contrary to the popular view. 

We studied the time evolution of the density for a variety of system sizes and initial densities 
(always taking equal concentration of A and B particles). We considered systems up to size 4 x 10^ 
for up to 5 X 10^ time steps. Some example runs are shown in figure |l[ We see that qualitatively the 
asymptotic behavior seems to be independent of the initial density, and that there appears to be a 
power law decay of the density with exponent approximately equal to —1/3. 

To get a better view of the slope of these curves we examined the behavior of p{t)/p{t/2) as a 
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Figure 2: Semilog plot of \og2{p{t) / p{t /2)) vs. time for system of size 4 x 10^. Initial densities are 
0.9, 0.5 and 0.2. 



function of time. Specifically, \og2[p{t)/ p{t/2)\ should give the slope of the line {i.e. the exponent of 
the power law behavior). Data for the same three runs as in figure H are plotted in figure |2[ This 
confirms that the decay behaves similarly for the different initial densities, although the noise is now 
much more apparent. To reduce the effect of this, we averaged the slopes illustrated in figure ^ over 
11 different initial densities. The resulting data is reproduced in figure ^. The errorbars represent 
the standard deviation of the different slopes. 

The errorbars at early times are primarily systematic in origin. They represent the fact that 
the decay rates are initially different for different initial densities — the particle motion is diffusive 
in nature until the particles have had a chance to interact with each other, so that we see t~^'^ 
behavior at early times for low density. On the other hand the errors at large time are primarily 
random errors; fluctuations in the initial density are significantly amplified when the density gets 
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Figure 3: Semilog plot of average slope vs. time for system of size 4 x 10^ 



small. 

Extrapolating the data of figure ^ to infinite time is problematic. It is not clear whether the 
effective exponent will continue to decrease; all one can say with confidence is that it is less than 
—0.31. There are arguments, however, for believing that the answer should be —1/3. 

The exponent —1/4 for the diffusive case arises because the relevant diffusive length scale at time 
t is t^". On the other hand, in the noisy Burgers equation (^) and thus in the one dimensional 
ASEP, there is super diffusive behavior, so that fluctuations spread more quickly and the length scale 
is larger. Thus the relevant length scale is t'^'^ [10, 11 1. Thus one might expect the concentration to 
go as 



c(t)~LVVL~ t2/3 ''yt2/3^^ 
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(5) 



certainly consistent with what is observed. Of course this argument is neglects the fact that only 
isolated clusters of one species undergo asymmetric simple exclusion dynamics, but one expects that 



clusters each species do spend most of their time reasonably isolated from each other ||5|, p, so that 
this analysis should be applicable. 
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